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Рассматривается характеристический полином, коэффициенты которого заданы числовыми интервалами. На основе многопа- 
раметрического интервального расширения метода корневого годографа определяются условия принадлежности образов вер- 
шин многогранника коэффициентов границам корневых областей. Разрабатывается методика нахождения проверочных вер- 
шинных полиномов для анализа минимальной степени устойчивости и максимальной степени колебательности системы с ин- 
тервальными параметрами. 


Введение 

В настоящее время в теории автоматического 
управления активно развивается робастное напра- 
вление. Одной из актуальных задач робастного 
управления является оценка качества автоматиче- 
ских систем, в которых некоторые физические па- 
раметры не определены точно или могут изменять- 
ся по заранее неизвестным законам в определен- 
ных пределах. Согласно [1-3], возможный путь ее 
решения состоит в использовании характеристиче- 
ских полиномов с интервальными коэффициента- 
ми, которые находятся по правилам интервальной 
арифметики. Для анализа системы на основе ин- 
тервальных характеристических полиномов (ИХП) 
главным образом применяются алгебраический и 
частотный подходы и в меньшей мере - корневой 
метод. Однако известно [4] , что робастное расши- 
рение корневого метода открывает новые возмож- 
ности для исследования качества систем с интер- 
вальными параметрами. 

Выбирая для решения поставленной задачи 
корневой подход, заметим, что поскольку интер- 
вальные коэффициенты имеют фиксированные 
пределы изменения, то корни ИХП мигрируют по 
комплексной плоскости, оставаясь локализован- 
ными в некоторых замкнутых областях. Если эти 
области находятся в левой полуплоскости, то по их 
расположению можно определить диапазоны воз- 
можного изменения корневых показателей каче- 
ства системы - степени устойчивости и степени 
колебательности. При этом особый интерес пред- 
ставляют их предельные значения - минимальная 
степень устойчивости и максимальная степень ко- 


лебательности. Именно они характеризуют наих- 
удшие режимы работы системы при изменении ко- 
эффициентов ее характеристического полинома. 

1. Постановка задачи 

Пусть ИХП имеет вид 

0(х) = а п х п +а п _у~' +--- + а 0 , ( 1 ) 

где коэффициенты я„ /=0,я могут изменяться в из- 
вестных пределах. Образуемый этими коэффици- 
ентами многогранник Р т является прямоуголь ным 
гиперпараллелепипедом с вершинами Ѵ р д=1,2”. 
Пусть каждой точке Р т т=п+1, соответствует 
устойчивый характеристический полином с по- 
стоянными коэффициентами и, следовательно, 
фиксированные корни в левой полуплоскости 
комплексной плоскости, определяющие степень 
устойчивости и степень колебательности системы. 
Необходимо оценить минимальную степень устой- 
чивости и максимальную степень колебательности 
для заданного многогранника Р т . 

Решить данную задачу возможно путем отобра- 
жения Р т на плоскость корней и нахождения обла- 
стей локализации корней ПХП. В соответствии с 
известной реберной теоремой [5] эти области огра- 
ничиваются образами ребер Р т . Поэтому одним из 
способов решения поставленной задачи может 
быть реберный анализ - отображение ребер Р т . Не- 
достаток такого подхода заключается в значитель- 
ной трудоемкости отображения всех (т- 2'"~ 1 ) ребер. 

Сократить число отображений возможно за 
счет применения предложенной в [6, 7] реберной 
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маршрутизации Р т . Она позволяет определить те 
ребра Р т которые непосредственно отображаются 
на границы корневых областей. 

Наряду с реберным анализом для оценки лока- 
лизации корней ИХП применяется также и вершин- 
ный анализ. Он основан на проверке расположения 
корней семейства полиномов, соответствующих 
вершинам Р т (вершинных полиномов). Так, соглас- 
но [1], для решения поставленной задачи необходи- 
мо найти образы всех 2" вершин Р т . Такая процеду- 
ра также является весьма трудоемкой. Следует заме- 
тить, что образы только некоторых вершин принад- 
лежат границам корневых областей (назовем такие 
вершины граничными). Поэтому для уменьшения 
вычислительных затрат представляет интерес задача 
определения граничных вершин. 

2. Графический способ определения 

граничных вершин 

Координ аты любой точки Р т относительно вер- 
шины Ѵ ѵ <7=1,2”, определяются выражениями 

а і = аі + До., / = 0, и; (2) 

(а — а‘ І )<Аа<(а -я в ), (3) 

где А а, - приращение /-го интервального коэффи- 
циента, о? - его значение в вершине Ѵ г Верхний 
предел коэффициента а, обычно обозначают сим- 
волом а,, а его нижний предел - символом а,. Со- 
отношение, связывающее точки Р т с корнями по- 
линома (1), может быть получено подстановкой в 
(1) выражения (2) 

0“(а) + ^Аа/ = 0, /=0Я (4) 

і 

где І)‘ ? (5) = ^а/У - вершинный полином. Введем 

І 

в рассмотрение ребра Р т , которые обозначим Я А 
где / - индекс интервального коэффициента, <? - 
индекс вершины, из которой по ребру изменяется 
а,. На основании (4) запишем уравнение отображе- 
ния К? на плоскость корней 

ОЧ.?) + АаД=(). (5) 

Согласно теории корневого годографа уравне- 
ние (5) может быть получено из реберной переда- 
точной функции 

Ж«(Д а,,5) = ^-. (6) 

' ' В«(з) 

Анализируя (5) и (6), заметим, что при измене- 
нии А а, в интервале (3) корни движутся от полюсов 
функции (6), соответствующих одному концу К?, к 
корням (5) на другом конце Д". При этом они обра- 
зуют фрагментарные ветви корневого годографа. 
Назовем их реберными ветвями (обозначим ЯА'% а 
их начала и концы - корневыми узлами ( Щ. 

Исходя из уравнения фаз корневого годографа 
[8], угол выхода ЯА‘‘ из комплексного корневого уз- 
ла ІІ Ч при увеличении а, находится по формуле 


Ѳ*=7Г-2Х+/Ѳ 0 , (7) 

а при уменьшении а, 

Ѳ ?=-ІЛ + /Ѳ о> (8) 

где Ѳ Е и Ѳ 0 - углы между вещественной осью и век- 
торами, направленными из ІІ д соответственно к 
#-ому полюсу и к / нулям передаточной функции 
(6) с координатами (О'фО). 

В [6] на основании (7) и (8) получено следующее 
условие отображения вершины Р т на границу кор- 
невой области: Д будет принадлежать границе обла- 
сти локализации корней, если разница между мак- 
симальным и минимальным углами Ѳ? меньше 180°. 

Заметим, что согласно (7) и (8) у всех углов Ѳ? 
есть общая составляющая |)Ѳ 4 , которая не влияет 
на выполнение условия граничного положения (Д 
Поэтому перейдем от рассмотрения углов Ѳ? к рас- 
смотрению углов Ф/, рассчитываемых при увеличе- 
нии а, по формуле 

Ф «=* + /©„, (9) 

а при уменьшении а, - по формуле 

Ф ? = /Ѳ 0 . (Ю) 

Введем в рассмотрение полярную систему ко- 
ординат с векторами единичного радиуса Д, выхо- 
дящими из начала координат под углами Ф?. Назо- 
вем векторы Е, коэффициентными векторами. 
Если задать значение Ѳ 0 , то согласно (9) и (10) угол 
вектора для каждого коэффициента легко опреде- 
ляется без вычисления корней соответствующего 
вершинного полинома. Очевидно, что по взаимно- 
му расположению Д можно судить о характере узла 
ІІ Г Он будет принадлежать границе корневой обла- 
сти, если векторыД всех интервальных коэффици- 
ентов будут располагаться в угле, меньшем 180°. 

Наиболее просто и наглядно проверить выпол- 
нение этого условия можно графическим способом, 
строяД на круговой диаграмме. На рис. 1 на приме- 
ре ИХП с четырьмя интервальными коэффициен- 
тами представлено возможное расположение Д для 
случая граничной вершины. Сплошной линией на 
рисунке показаны векторы при максимальном зна- 
чении интервального коэффициента в узле (Д а 
прерывистой - при его минимальном значении. 



Рис. 1 . Коэффициентные векторы граничной вершины 
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Заметим, что изменение предела коэффициента 
с минимального на максимальный или наоборот 
позволяет поворачивать соответствующий вектор 
на 180°. Это дает возможность располагать векторы 
желаемым образом (в соответствии с условием гра- 
ничности корневого узла). Получающиеся пределы 
коэффициентов ИХП будут определять координа- 
ты граничной вершины Р т . 

На основании анализа возможного расположе- 
ния Е і можно сделать следующий вывод: если для 
каждого из т коэффициентов ИХП построить на 
круговой диаграмме в соответствии с (9) и (10) два 
противоположных вектора (рис. 2), то любые т по- 
следовательно расположенных векторов будут ле- 
жать в угле, меньшем л, и, следовательно, опреде- 
лять координаты граничной вершины Р т . Задавая 
направление вращения и выбирая в качестве на- 
чального каждый из 2т векторов, можно получить 
набор из 2т возможных граничных вершин Р т . 



Рис. 2. Диаграмма расположения коэффициентных векторов 

3. Анализ влияния особых лучей 
на граничные вершины 

Установлено [6] , что узлы Ѵ я , граничные в одной 
области комплексной плоскости, могут становить- 
ся внутренними в другой области. Это свойство 
определяется характером границ областей локали- 
зации корней ИХП [7]. Реберные ветви КО," могут 
входить в состав границ как полностью, так и ча- 
стично, пересекаясь в определенных точках. Для 
анализа данного свойства рассмотрим две гранич- 
ные реберные ветви с общим корневым узлом, яв- 
ляющиеся образами двух ребер грани Р„. Обозна- 
чим грань через Су, где сі и /- индексы изменяю- 
щихся на грани интервальных коэффициентов. 
Согласно [7], если образы ребер грани 0^ пересека- 
ются в точке Ф=а+іР, то координаты У являются 
решением системы уравнений 

' МСО'ПтДУ/ЬКеД/уОітДУ)^ =0; 

< Яе((У/)1т[Х«^*)^+Х«^*) Р ]- 

к р 

. -Іт((*У)МІ >*Ѵ)‘ + Х>,СО'] =0. (П) 

к р 

где к - индекс интервальных коэффициентов, при- 
нимающих граничные значения в вершине ц. 


Система уравнений (11) записана для общего 
случая ИХП, когда в нем есть и постоянные коэф- 
фициенты с индексом р. На основе формулы Муав- 
ра $ ь =\У(со$(Ь(р)+ріп(Ь(р)) преобразуем (11) к виду 

" I + \сІ + / 

я (со&(сіф)&т(/ф) -со&(/ф)$'т(с1ф)) =0; 

< Л а к У | /+А 5ІП((* -/» + 

к 

+Х а 1 1 5 * | 7 +Р 8 ІП ((Я - /)<Р) = 0- (12) 


В (12) координатами У являются ф и |з’|. Так как 
|Т|^0, первое уравнение системы (12) можно упро- 
стить до вида 

8Іп((/-йО<р)=0. (13) 


Решением (13) при (/-Д)ф0 и учете периодично- 
сти синусоидальной функции является выражение 


Ф = 


ЛІ 


/= 0 , 1 , 2 ,... 


(14) 


Из ( 14) следует, что может быть несколько точек 
У пересечения реберных ветвей грани Оу. Получае- 
мые из (14) углы ср задают особые лучи, на которых 
У могут находиться. Если обозначить \ф—с1\ =?„ тогда 
для любой грани Р т углы особых лучей определя- 
ются по формуле 

<Р- =— , г = 1,2,3,..., и. (15) 

2 


Таким образом, при определении набора гра- 
ничных вершин Р т необходимо знать углы ф 7 всех 
особых лучей для заданного ИХП. Очевидно, что 
для анализа корневых показателей качества устой- 
чивой интервальной системы следует выбрать ф г из 
диапазона 90°<^,<180° и определить образуемые 
ими во втором квадранте секторы. Далее для каж- 
дого из этих секторов необходимо построить круго- 
вую диаграмму коэффициентных векторов и опре- 
делить соответствующий набор граничных вер- 
шин. Полное множество граничных вершин Р т яв- 
ляется объединением множеств вершин, получен- 
ных для различных секторов. 


4. Методика определения граничных вершин 

Как было указано в постановке задачи, для ана- 
лиза минимальной степени устойчивости и макси- 
мальной степени колебательности системы с за- 
данным ИХП достаточно определить корни поли- 
номов только в граничных вершинах Р т . На осно- 
вании проведенного выше анализа разработана 
следующая процедура формирования набора гра- 
ничных вершин для левой полуплоскости корней. 

1. Проверить робастную устойчивость ИХП (на- 
пример, с помощью полиномов Харитонова [1]). 

2. На основании (15) найти углы особых лучей, 
принадлежащие диапазону 90°<^,<180° и опре- 
делить секторы, на которые лучи разделяют 
второй квадрант плоскости корней. 
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3. Из каждого сектора произвольно задать по од- 
ному значению угла Ѳ 0 . 

4. Для минимального Ѳ 0 на круговой диаграмме 
для всех интервальных коэффициентов ИХП с 
обоими пределами построить векторы под угла- 
ми, рассматриваемыми по формулам (9) и ( 10). 

5. Начиная с любого коэффициентного вектора, 
выбрать т векторов, расположенных последо- 
вательно в положительном направлении (про- 
тив часовой стрелки) и определить координаты 
соответствующей векторам граничной верши- 
ны Р т . 

6. Последовательно изменяя начальные векторы, 
повторить процедуру п. 4 и получить координа- 
ты 2т граничных вершин Р т . 

7. Повторить действия пп. 3-5 для возрастающих 
значений выбранных в п. 2 углов Ѳ 0 . При каж- 
дом последующем угле Ѳ 0 из получающихся на- 
боров граничных вершин выбирать новые вер- 
шины и добавлять их к найденным ранее. 

5. Пример 

Пусть система имеет характеристический поли- 
ном Д(5)=а 3 і 3 +й 2 5 2 +й 1 5+й 0 , где й 3 =[ 1 ; 1,2] , 
й 2 =[ 12; 12,2] , й,=[22;23,2], й 0 =[20;25,5]. Требуется 
определить минимальную степень устойчивости и 
максимальную степень колебательности интер- 
вальной системы. 

Для решения этой задачи вначале было устано- 
влено, что заданный ИХП робастно устойчив. Да- 
лее согласно разработанной методике найден один 
особый луч с углом (р 3 =120°. Он разделяет второй 
квадрант на два сектора, из которых выбраны два 
угла: Ѳ 0 =150° и Ѳ 0 =1ОО°. Для первого угла постро- 
ена круговая диаграмма коэффициентных векторов 
и определен следующий набор граничных вершин 
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й, , й, . й| . й„ , @ 2 , й 2 , й 3 , Йо , й 3 , й 2 , й 3 , йд , 


а 3 ,а 2 ,а ѵ а 0 \ а 3 ,а 2 ,а х ,а 0 ; а 3 ,а 2 ,а х ,а 0 ; а 3 ,а 2 ,а х ,а 0 . 

При Ѳ 0 =ЮО° по второй диаграмме векторов до- 
полнительно найдены еще две граничные вершины 


Для каждого полученного набора коэффициен- 
тов определены корни соответствующих полино- 
мов. В результате анализа корней установлено, что 
рассматриваемая система имеет минимальную сте- 
пень устойчивости а тіп =0,95 в вершине а ь а 2 ,а х ,а 0 и 
максимальную степень колебательности ц шах = 1,37 
в вершине а ъ а ъ а ь а п . 

Заключение 

Для оценки корневых показателей качества си- 
стемы с интервальным характеристическим поли- 
номом разработана методика определения гранич- 
ных вершин многогранника интервальных коэф- 
фициентов. Образы этих вершин принадлежат гра- 
ницам областей локализации корней интервально- 
го полинома и определяют минимальную степень 
устойчивости и максимальную степень колебатель- 
ности системы. Применение разработанной мето- 
дики позволяет сократить число проверяемых вер- 
шин и, следовательно, упростить процедуру анали- 
за. В основу методики выбора граничных вершин 
положены фазовые соотношения известного мето- 
да корневого годографа. Предложенный графиче- 
ский способ определения координат граничных 
вершин отличается простотой и наглядностью. 
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